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We present a numerical study of a two dimensional model of the Wess-Zumino type. We formu- 
late this model on a sphere, where the fields are expanded in spherical harmonics. The sphere 
becomes fuzzy by a truncation in the angular momenta. This leads to a finite set of degrees of 
freedom without explicitly breaking the space symmetries. The corresponding field theory is ex- 
pressed in terms of a matrix model, which can be simulated. We present first numerical results 
for the phase structure of a variant of this model on a fuzzy sphere. The prospect to restore exact 
supersymmetry in certain limits is under investigation. 
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SUSY an a fuzzy sphere 



1. The Di Vecchia-Ferrara model 

Back in 1977 Di Vecchia and Ferrara introduced an elegant supersymmetric (SUSY) model of 
the Wess-Zumino type in two dimensions [|l|]. On a Euclidean plane its action reads 



S[(j),xif) = \ jd 2 x [d^d^ + yy v d v y + [V'm 2 + WV"(0) 



(1.1) 



where is a scalar field, and \j/ is a 2-component Majorana spinor field (with y = \j/ T C, where 
C = ic% is the charge conjugation operator). 

We fix the boson-fermion interaction, as well as the bosonic potential, by the choice 



v <« = Vl^ 3 + 7* 2 <12) 

This action is invariant (up to a total divergence) under the SUSY transformation 



50 = ey 

Sy=[y v d v $-V'(<!))}e , 8y = s[-y v d v <p - V'(<j>)] 



(1.3) 



where £ is a constant spinor field. 

2. The fuzzy sphere 
2.1 Geometry 

The coordinate operators x, on a fuzzy sphere of radius R have to solve the equation |Q] 

x 2 l +x 2 2 +xl=R 2 \ . (2.1) 

This can be achieved by identifying the position operators Jc; with the (re-scaled) angular momen- 
tum operators L,- , 

R li, i = 1,...,3 , M , (2.2) 



where I is the spin of the irreducible representation. For finite £ we obtain N x N matrices, with 
N = £ + 1 . Then the coordinate operators on the fuzzy sphere are non-commutative, 

£ijA ■ (2.3) 



In the limit I —> oo commutativity is restored, and the sphere is not fuzzy anymore. 
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2.2 Fields 

Some field, for instance our scalar field <f>(x), on the sphere S 2 can be expanded in term of 
spherical harmonics Y(. m {x) , 

oo e 

<H*) = L L C im Y £m(x) ■ (2-4) 

(=0 m=-£ 

In accordance with the above treatment of the coordinates, we are going to represent also the 
fields as matrices, so that we end up with a matrix model. Without limiting the matrix size N, a 
Hermitian matrix can be expanded in the polarisation tensors Yf m [||] , in analogy to eq. (2.4), 



= L L c e>nYe m • (2.5) 

£=0 m=-£ 

The polarisation tensors play a role analogous to the spherical harmonics. In particular they obey 
L 2 Y( m = £(£+ l)% m , where L 2 is the angular momentum squared (cf. Section 3). Moreover we have 
the adequate normalisation and parity behaviour, ^Tr(Yj m ,Y£ m ) = 8e>t8 m r m , Yj m = {—\) m Yt_ m . 

We now introduce a cutoff £ max = N — 1 , and the remaining coefficients for the field <p can 
be embedded into Hermitian N x N matrices. Thus we arrive at a finite set of degrees of freedom, 
without any explicit breaking of the space symmetries. Therefore this regularisation is attractive for 
SUSY models, where the lattice formulation is notoriously troublesome Moreover the fuzzy 
sphere regularisation is not plagued by the fermion doubling problem [j5j]. 

Theoretical aspects of SUSY on a fuzzy sphere are discussed in Refs. ^ and reviewed in 
Ref. [0]. The feasibility of numerical simulations in this regularisation has been tested for the A0 4 
model in d = 2 [j|] and d = 3 [|h. For an approach to simulate SUSY gauge theory on the fuzzy 



sphere, see Ref. [10] 



3. The Di Vecchia-Ferrara model on a fuzzy sphere 

To be explicit, we transfer the Di Vecchia-Ferrara model from the Euclidean plane to a fuzzy 
sphere by means of the following substitutions: 

i r . , , i 



3 

Mix) -> U = 'AUM , -d 2 Hx) -> (3-D 



where <p is a Hermitian N x N matrix. We are ultimately interested in the limits N, R —>■ °°. 

In practice the left-handed and right-handed applications of the operators can be implemented 
best by storing the matrix configurations (j> as vectors. In this setting the Dirac operator ( JO) ) takes 
the form of a 2N x 2N 2 matrix. We symmetrise the fermionic potential as 

Kym(^) = \(v"($)®lL N +Tl-N®V"($f) , (3.2) 

which leads to the Dirac operator 

D= (d 3 + ±1 N + V» m ^ \ 
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phase 


<0,?> 


(I0I 2 ) 


disordered 


rafj 


«0 


uniform ordered 


>0 


>0 


non-uniform ordered 




>0 



Table 1: The phases that we observed, along with the respective magnitudes of the order parameters. 



with d± = d\ ±ic?2, and with the representation (3.1) of the differential operators. As usual in 
Euclidean space we deal with an anti-Hermitian kinetic part. This also includes the term i/R, 
which emerges from the curvature effect on the spin connection. 

Integrating out the fermionic variables yields the Pfaffian Pf [(CD) asym ], where the subscript 
means the anti-symmetric part. As a first approach to explore this type of model, we replace the 
Pfaffian by VdetD, 1 so we arrive at the matrix model action 

m = ^{\t * ]] + t* 2+ m \fr2* 3 + t\ ^ 4 ) " ln v ^ • (3 - 4) 

4. Order parameters 

As in Refs. [Q, ^] we are going to explore the phase diagram by considering order parameters, 
which are constructed from the coefficients C£ m in the expansion ( |2.5[ ). They can be extracted from 
the relation 

Cem = —Tr[(Y em )^}. (4.1) 

In particular we focus on the quantities 

1 An 
$ ■= £ hm\ 2 and |0| 2 := £> 2 = — Tr [f] . (4.2) 

m=-i I ly 

Based on the magnitudes of the expectation values and (\(j)\ 2 ) we distinguish three phases as 
specified in Table [j], in close analogy to the A0 4 model on a non-commutative flat space [O]: 

• In the disordered phase (tyf) holds for all i. The angular mode decomposition does not 
detect any contribution that could indicate spontaneous symmetry breaking. 

• The uniform ordered phase is characterised by (|0| 2 ) ~ (0 ( 2 ) 3> 0, i.e. the zero mode con- 
tributes significantly, whereas higher modes are suppressed. This phase corresponds to the 
spontaneous magnetisation in an Ising-type spin model. 

• In the non-uniform ordered phase a non-zero mode condenses. This leads to the relations 
(I0I 2 ) « (l<H 2 -0o> »°- Inthis 

case the rotation symmetry is spontaneously broken. That 
phase is specific to the fuzzy sphere; it does not occur in commutative spaces. 

We add that the fuzzy sphere formulation is non-local, as the non-commutativity of the coordi- 
nates suggests. Therefore the Mermin-Wagner theorem does not prohibit the non-uniform ordered 
phase. 

'The impact of this substitution remains to be investigated in detail. Further comments are added in Section 5. 
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5. Numerical results for the phase diagram 



In general the determinant detD (where D is given in eq. ( p.3| )) is complex. In the final limit 
N, R — ► °° it is real positive, however, hence the complex phase represents an artifact of the fuzzy 



regularisation (apart from the substitution in eq. (3.4)). We therefore modify the regularisation at 



this point by using |detD| already at finite N and R. With this modification the action (3.4) defines 
a Boltzmann weight which enables Monte Carlo simulations. We performed such simulations with 
the Metropolis algorithm: in each step, a conjugate pair of matrix elements is updated, and |detD| 
is computed explicitly. Throughout our simulations we fixed the radius of the sphere to R = 1. 

So far we simulated at N = 4, 6 and 8 (which are numerically handled by matrices of size 
32, 72 and 128) and we explored the phase diagram in the (m,A) plane. We also measured the 
phase of the determinant. It turns out to be quite stable, which is favourable for the modified 
regularisation (in the extreme case of an invariant phase the modification is redundant). Examples 
for this property are illustrated in Fig. |l[ 



N = 6, 1=4, m =4 
N = 8, % = 1, m = -6 



200 400 600 800 1000 1200 1400 
configuration number 

Figure 1: Two examples for histories of the phase of the fermion determinant. The parameters in these 
simulations are N = 6, A = m = 4 and N = 8, X — 1 , m = —6. 

Next we show in Figs. || and || the expectation values of the order parameters (|0| 2 ), (Qq) 
and (|0 1 2 — 0q) (as discussed in Section 4) for N = 6 and N = 8, at X = 1, 5 and 9. We actually 
investigated a larger range of the mass parameter m , but we show here the interval of interest in 
view of the phase diagram. In all cases, large values of \m\ lead to the disordered phase. When 
\m\ decreases below a critical value 2 (which is similar but not identical for both signs of m) we 
enter the phase of uniform order. In the vicinity of m = we also observe the phase of non-uniform 
order to set in; the corresponding conditions are given in Table [jj 

For N = 6 and 8 we probed A = 1 , . . . , 10 , which gives rise to the phase diagram in Fig. |]. 



6. Conclusions 

We explored a new way to simulate a two dimensional model of the Wess-Zumino type. The 
model is wrapped on a sphere and the fields are expanded in spherical harmonics. A truncation 
in the angular momentum renders the sphere fuzzy, and the corresponding field coefficients build 
a finite set of degrees of freedom to be used in numerical simulations. In this first approach we 
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Figure 2: The order parameters for N = 6 and A = 1,5,9 in the range m = — 8 . . . 8. We see a double 
peak structure in all cases, and the peak at m rs —2 is enhanced for increasing A. A strong coupling A also 
implies a gradually enlarged phase of non-uniform order. 



>l 2 






< lif > ' ' ' ' 


1 = 5 : 


< 4> 2 > l - * - j 




<W Z -t Z > l-B-l 








— * — i — * -* , — i — « — , 




Figure 3: The order parameters for N = 8 a«<i A = 1, 5, 9 in the range m = —8 ... 8. The phase transitions 
show a behaviour similar to the results for N = 6 ( Fig. 2 ), although the values of the order parameters differ. 
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Figure 4: The phase diagram at N = 6 and 8, which is identified from measurements as shown in Figs. 2 
and 3 for A = 1,2, ...,10. 

simplified the Pfaffian to V detD. Thus we studied a SUSY inspired system of interacting scalars 
and Majorana fermions on a fuzzy sphere. 

In the final limit of infinite angular momentum cutoff ,/V and radius R, the determinant detD is 
real positive. This does not hold at finite N and R, so we modify the regularisation by employing the 
modulus of the fermion determinant already on the regularised level. We expect this formulation 
to lead to the same limit. This expectation is supported by the observation that the fluctuations of 
the phase arg(detD) are small. 
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With this method we simulated the system at N = 4, 6 and 8 and R = \. The basic properties 
are similar in all cases; in particular large \m\ >2 always leads to a disordered phase. So it is con- 
ceivable that we are already peeping at aspects of the large N limit. However, the final stabilisation 
of the phase diagram at large N may involve a re-scaling of the axes. 

The ordered non-uniform phase emerges as a consequence of the non-commutativity of the 
coordinates, which we use on the regularised level. We observed that phase around m pz 0, and it 
ought to evaporate as we proceed to larger N. That feature is left for further investigation. 

Although this project is still on-going, the preliminary results are encouraging regarding the 
hope to find a way to formulate and explore SUSY inspired models beyond perturbation theory. 
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